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Abstract 

An  electric  current  source  comprising  point  horizontal  electric  dipoles  is  proposed.  An 
approximate  description  of  the  magnetic  field  from  such  a  source  is  given.  The 
environments  con^ered  are  two-layer  (2  semi-infinite  media,  typically  the  air  and  the  sea) 
and  three-layer  (2  semi-infinite  layers  separated  hy  a  conducting  slab,  typically  the  air,  sea 
and  sea-bed). 
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Approximate  Expressions  for  the 
Magnetic  Field  from  an  Electric  Current 
Source  in  a  Layered  Conducting 
Medium 


1.  Introduction 

Reference  1  describes  an  electric  source  mathematical  model  in  which  a  set  of 
electric  source  parameters  can  be  deduced  from  a  measured  electric  field  matrix. 
The  enviroiunent  is  modelled  by  two  semi-infinite  media  (the  air  and  fite  sea¬ 
bed,  the  former  having  an  electrical  coixiuctivity  of  zero)  separated  by  a 
conducting  slab  of  finite  thickness  and  non-trivial  conductivity  (foe  sea).  The 
source  is  modelled  as  either  a  line  of  point  current  sources  and  sinks,  or  by  a 
continuous  distribution  of  these,  the  density  distribution  function  being 
described  by  a  polyrtontial.  In  either  case,  the  total  current  must  be  zero.  In 
the  discrete  case  the  amplitudes  of  the  currents  must  sum  to  zero,  whilst  in  foe 
continuous  case  the  integral  of  the  distribution  function  must  be  identically 
zero. 

The  effectiveness  of  the  point  source  rrvxiel  to  describe  foe  electric  fields  has 
been  well  showiu  However,  there  is  considerable  interest  in  determining  foe 
magnetic  fields  from  such  a  source. 

The  electric  field  is  determined  by  (rdatively)  simple  image  theory. 

However,  aj^lication  to  the  magnetic  field  case  is  not  as  straightforward  smce 
the  currents  (which  produce  foe  field)  are  distributed  through  foe  volume  of  the 
conducting  media. 

In  this  report,  the  simpler  case  of  just  two  semi-fofinite  media,  with  a 
common  bouivlary  (typically  foe  air  and  an  infinitely  deep  sea)  is  considered 
first  Then  foe  more  general  three-layer  case  is  cortsidered.  In  both  cases,  the 
aim  is  to  produce  eiqnessions  for  the  magnetic  field  from  a  point  dipole  source 
rather  fow  current  souroe/sinks.  Smce  foe  latter  must  sum  to  zero  to  be 
jfoysically  realizable,  any  such  source  of  n  current  source/sinks  can  be 
represented  by  (n-1)  dipoles. 


2.  Two-Layered  Medium 


Consider  the  geometry  of  Figure  1,  in  whidt  layer  1  is  infinitely  deep  (d  =  «>),  or 
altemativdy  Oj  » In  pr^cal  terms,  fiiis  is  equivalent  to  very  deep  sea¬ 
water,  where  the  source  a^  sensor  are  not  near  the  bottom.  For  a  point 
horizontal  dipole  in  the  x-axis,  of  amplitude  P,  at  (Xq,  y^.  Zq).  refdcnce  2  gives 
the  expected  magrtetic  fields  in  layer  0  and  layer  1  at  die  point  (x,  y,  z).  These 
are  given  in  cylindrical  coordmates  as  follows  (for  the  case  where  Oq  =  0): 


It  is  of  particular  interest  that  B  is  indepeitdent  of  either  conductivity  (o^y  O]), 
except,  of  course,  diat  the  conductivity  of  the  layer  containing  the  source  will 
determine  the  current  flow  (hence  P),  at  least  for  a  corrstant  voltage  source  (as 
would  be  the  case  for  a  galvanic  source  such  as  the  potential  differerKe  between 
a  ship's  propdlor  and  hull). 

1b  convert  die  dwve  to  Cartesian  coordinates,  the  following  appl)r: 
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Bx  =  6,  cos  0  -  B^  sin  ^ 
By  =  B,  sin  ^  +  B^  cos  ^ 


and 


sin^  =  (y  -  Vo),  cos  ♦  =  (x  -  a[<P,  =  (x  -  Xq)*  +  (y  -  yo)^ 


Figure  1:  Geometry  Definition 


! 
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3.  Three-Layered  Medium 

Now  consider  the  full  geometry  of  Figure  1.  In  this,  the  three  layers  have 
coiKluctivities,  Oq,  o,  and  c^.  Typically,  for  medium  0  being  air,  Oq  =  0.  A 
source.  So  is  located  at  (x^  yp,  Zq)  within  medium  1,  whose  thickness  is  did  >  0). 

In  the  first  instance,  this  source  is  taken  as  a  point  source  or  sink  of  current,  I, 
aiKi  is  considered  in  isolation.  RefereiKe  3  gives  the  electric  scalar  potential,  V, 
in  the  three  media  as  follows: 


I 


where 

r2  =  {x-xoi^  +  iy-y^ 

Qj  =  (aj  -  Ofl)  /  (o,  +  Oq),  (32  =  (o,  -  aj)  /  (o,  +  Oj) 

Following  reference  4,  the  magnetic  field  from  this  source  is  determined  by 
invoking  Ampere's  law 

V  A  B  =  po 

Thus,  by  Stokes'  theorem: 

J  V  A  B.dA  -  jB.dl  -  P50  jE.dA 


This  problem  can  be  considered  as  one  with  cylindrical  symmetry,  with  the 
z-axis  being  the  axis  of  symmetry.  Consider,  then,  the  geometry  of  Hgure  2, 
where  the  cable  carrying  current  to  this  ’monopole"  source  is  explicitly  iircluded 
(the  field  from  this  is  subtracted  later). 


i.e.  Ej  is  independent  of  ^ 


J"  E .  dA  m  2n  j  rdr 


.Inr  «  PoO„2it  j E^rdr 
0 


where  n  refers  to  the  medium  containing  the  field  point  (n  =  0, 1,  2). 


Now,  Ej  =  -  3V  /  3z,  and  V  is  known.  In  general. 


EE. 


where  EJ  is  of  the  form: 


rdr 

[r^  +  (z 


and  the  [j]  refer  to  the  various  source  and  image  terms. 


f  E  rdr  -  __L 
J  h  4iio 


•••  \  ■ 


I 

4no 

1*0^ 

4«r 


(Z-Zj) 


(z  -Zj) 


sgn(z-z.)  -  — - - '  . 

[r^  +(z  -z.)^]’^* 


(z  -Zj) 


12 


Let  Oq  =  0  Qi  =  1  and  Q2  can  be  written  as  Q.  For  medium  0  (air) 

z>0,-d<ZQ<0. 

Since  a  o^  then  =  0. 

However,  the  field  horn  the  supply  cable  must  be  subtracted  from  this  6^. 

For  this  medium,  we  choose  the  cable  to  run  from  -  oo  to  Zq.  The  field  from  the 
cable  in  this  instance  is  given  by: 


B. 


1^^ 

4nr 


(z-Zq) 


\  ■  — 
*0  4nr 


1 


(z  -z.) 


(r^*(z-Zg)^y'^ 


Note  that  the  field  in  air  is  quite  indepotdent  of  water  depth  and  conductivity. 
Consider  medium  2  (sea-bed)  next,  where  z<-d, -  d<ZQ<0.  Here, 


sgn(z±Zo)  =  -1 

sgn(z-2kd±z^)  =  -1,  Vi:>l 


Hence,  implicitly  including  the  feed  cable,  B^  is  given  by: 


-Pof 

B.  =_^(1-Q) 
’>  4jtr 


1  + 


(z-Zq) 

(r^  +(z -Zaf)^'^ 


1  + 


(z  +Zn) 


V 


(r^+(z+Zo)V* 


E  Q" 


k-1 


1  + 


(z -2kd  -Zp) 
{r^*(z-lkd-z\Y?'^ 


-E  0“ 

k>l 


1  ♦ 


(z  -2faf  +Zq) 
(r^+(z-2jU+Zp)^)’'^ 


In  this  medium,  we  choose  the  cable  to  run  from  +  00  to  2^  for  which  the  field 
is  given  by: 


B. 


4nr 


(1  + 


(2-Zn) 


(r^+{z-2o)2)'/2 


Thus,  reduces  to; 


-M 

4nr 


1  - 


(Z-Zg) 

(r^  +  (z 


+  (1-Q)  E  0“ 

k<0 


( 

(.z-2kd  -Zg) 

(z-2JW+Zg)2)^/2  1 

(r^  +  (z  -Hcd  -z»)^)’^^ 
V 

(r^+(z-2JW+z.)2)’/2 

jj 

In  medium  1  (sea),  -  d<z<0,  -  d<Zg<0. 
Here, 


sgn  (z  +2kd±z„)  =  +1 

sgn  (z-2kd±za)  =  -1  >  ^ 

sgn  (z+Zg)  =  -1 


The  regions  z  >  Zq  and  z  <  Zq  must  be  considered  separately.  In  particular,  for 
z  >  Zj  the  feed  cable  is  assumed  to  run  from  -  <»  to  zq  and  sgn  (z-  Zg)  =  +  1. 
However,  for  z  <  Zg,  the  cable  runs  from  +  »  to  Zg  and  sgn  (z  -  Zg)  =  -  1.  In 
both  cases,  B^^  reduces  to; 


♦>  4nr 


(Z+Zg) 


*  (r-*(z*zg)^y'^ 

^  Q*(z-2fai-zg) 

M  (r^*(z-2kd-Zg)y^^  M 
^  Q‘(z*2W*Zg)  - 

w  {r^  *(z  *2M  *  w 


Q^(z-'2M  +Zg) 

(r^+(z-2W+Zg)^)’^^ 
Q*(z  +  2W-Zo) 

(r2+{z  +  2iW-Zg)2)’/2 


Since  the  total  current  must  sum  to  zero,  there  will  always  be  more  than  one 
source/sink  for  (his  soiuce  type.  Consequently,  it  is  convenient  to  represent 
the  source  as  a  line  of  point  (Upwles  as  in  reference  1.  The  magnetic  field  from 
one  such  point  dipole,  in  the  more  useful  Cartesian  coordinates,  is  determined 
as  follows. 


I 
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In  general,  we  may  say  that: 


■  -  B^  sin  <t).  By  »  B^  cos  (jt 


sin  <t)  -  (y  -Vq)  It,  cos  $  -  (x  -Xq)  It 


Let  there  be  monopoles  of  amplitude  -t-  /  at  (x  +  5x)  and  -  i  at  (x  -  Sx),  then  the 
field  from  such  a  dipole  is  given  by: 


where  X  =  x  -  Xq,  V  =  y  -  yo-  Spedfically, 

B  =  B  -  B 

y  ♦.  ♦. 


where 

rf  -  (X±8x)='  + 
B^^  -  B^(X±5x) 


Now,  contains  terms  of  the  form: 


At. 


Let  +  Zj^.  To  first  order. 
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1  1  4xbx 

“I  ■  “5  "  ~~T~ 
r*  r.  T 


I  1.1 


(X^-8x)  _  (X-Sx)  ^  28x  j_2  f 


(X+8x)  _  (X-8x)  ^  28x  ,  X  _  X 

r^rNzjV  '  r?(rSziV  '  [  '  l''J 


The  helds  may  now  be  determined.  However,  the  field  from  the  return  path 
must  be  included.  This  is  given  by  the  familiar  equation  for  a  point  electric 
dipole  source  in  a  homogeneous  medium. 

Specifically,  the  field  from  a  current  element  flowing  between  (X  +  8x)  and 
(X  -  8x)  is  given  by 


Po^  (X  +  8x)  _  {X-8x) 

^  1  (p2  +  (X  +  8x^)’''^  (p^+(X-8x)^)'''^ 


where  =  (y  -  yo)^  +  (z  -  Zq)^. 

9  is  in  the  (y  -  z)  plane  such  that  in  sin  0  =  (z  -  Zpl/p,  cos  0  =  (y  -  yo)/p. 
Expanding, 


4i^  r  hf 


where  P  =  2/8x,  the  dipole  moment  Thus, 
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-PoP<z-2o>  B  jioP(y-yo> 
iwF  4jiR^ 


Note  that  this  return  path  represents  the  only  source  of  vertical  field  and  is 
independent  of  depths  and  conductivities. 

Thus,  the  complete  dipole  fields  may  be  determined  as  follows: 

Medium  0 

Y 

B  2- 

*0  4ji  r  r  (r^  +  (2-2o)^)’'^  r^ 

y  ,  ^  f  ^ 

4jc  [  [f  J  J  (r^  +  (2 (r^ +(2 -2^)^)^^^  [ 

h  4jC  r3 


where  the  ±  should  be  read  as  "plus  and  minus",  implying  4  summations  in  this 
case. 


-\hp 


rr  , 


Xi 


f  V 
X 

v”y 


1-2 


1  +. 


(2 +20) 

(2+2o)  (2-2q) 

(r^+(z  +2q)^)*^^ 


7^ 


MqP  Y 

"JiT  "P 


Q‘'(z±2W±Zo) 

1  -2 

X 

T 

x^ 

(r^  +  (2  ±  2JW  ± 

/> 

(r^*iz±2M±z^f) 

Medium  2 


g  ^  X  y 

*j  4jt  r  r 


(2-2o) 


(r^  +(z  -Zq)^)^^^ 
-  'q*(2-2W±Zo) 

*  -Q>  r  — ! - Wj 

ft  (r^  +  (z-2JW±2o)^)’'^ 


(z-2o) 


(r^  +(2  -Zp)^)' 


7^  (r^*{z-2kd±z^f) 


B  .2^ 

4«: 


77 


1  -2 


X 

v”y 


1  - 


(z-Zo> 

(r^+Cz-Zo)^)’/" 


f 

w 

\ 

(  -] 

1 

1  - 

X 

/  > 

> 

B 


(2-2o) 

(r^  +  Cz-z')^)^'^ 


-  Q*(z-2W±2„) 

+  (1  -Q)  F  ^t-- _ - _ * 

M  {r^*{z-2kd±z^f)^'^ 

,  PqP  Y 
*2  4jt  r3 


f 

1  -2 

X 

r 

K  y 

(r^+(z-2JM±z.)2) 


The  above  expressions  can  be  readily  checked  by  confinning  the  continuity  of 
the  components  across  the  boundaries  at  z  =  0  and  z  =  -d.  Further  it  is  easy  to 
show  that  for  Q  -»  0  or  d  -»  ao,  the  expressions  reduce  to  those  given  earlier  for 
the  two  layer  case. 
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At  r  =  0  (directly  above  or  below  the  source)  the  expressions  as  given  are 
inappropriate.  However,  in  such  a  case  they  may  be  reformulated  as  follows: 


-  0  V  n  -  0, 1, 2 
B,  -  0  V  n  -  0, 1, 2 

*»|r^ 


Now, 


sgn  (z) 
sgn(z) 


The  following  are  obtained: 


/ 

1.1 

r 

2 

"z 

y 

1-^ 

r 

1  ^ 

z 

V  / 

L  to  first  order 


B. 


^-0  8k(z-Zo)^ 


1 


2-^0 


9 

I 


''  Y  ^ 


z  -z. 


/  J 


y. 


Ir-a 


ihp 

8tc 

t  Q‘ 

k«l 


(Z+Zg)'' 


-I 


.  V 


(z  +  2  fai  ±  Zp)^ 


1  - 


Z  +z„ 

X' 


V  °/ 
2 


2sgn(z  -Zq) 
(2-2o>^ 


(z  -2W±Zn)^ 


1-2 


4  (z  +  2iW±Zo)^  4  (z*2kd±z^f 
v2  7  v2  M 


4  {z-2kd±z^f  4  (z-2jM±Zo)2 


-l^oP 

• 

1  j 

/ 

1-31 

r  X 1 

■_  9 

Y 

8ic 

(Z-Zg)^ 

4 

r^J 

~  4 

+  (1-Q)5: 

3  Y* 


Q'‘ 


(z-2W±Zp)* 


9  X* 

4  iz-2kd±Zgf 


*(z-2kd±z^fj 


4.  Comparison  of  Two-Layer  and 
Three-Layer  Solutions 


The  above  expressions  for  the  two-layer  and  three-layer  case  are  readily 
progranuned.  IWo  sources  are  considered  here  as  examples. 

In  the  first,  a  horizontal  dipole  of  125  A.m  (a  current  source/sink  pair  of  50  A 
separated  by  25  m)  is  placed  at  3.35  m  depth  in  21  m  of  water  with  a  sea/sea- 
b^  Q-factor  of  0.8.  The  field  is  calculated  at  20  m  depth  for  the  two-layer  aiul 
three-layer  case  (directly  below  the  source)  and  the  results  are  shown  in 
Rgure  3.  In  this  geometry,  the  fields  are  entirely  athwartships. 

Qearly  the  use  of  two-layer  equations  would  produce  a  considerable  over¬ 
estimate  of  the  fields. 


-100  -60  -60  -40  -20  0  20  40  60  60  100 

Longitudinal  distance  Imetre) 


Ffgmv  3;  Bdm-xoater  magnetic  fidds,  simple  source. 


In  addition,  die  source  defined  in  Figure  3  of  Reference  1  (reproduced  here  as 
Figure  4)  is  also  used,  with  the  same  geometric  and  environmental  parameters. 
This  source  comprises  31  sources/sinks  (where  the  total  current  sums  to  zero) 
separated  by  2.67  m  and  has  a  fai^field  dipole  moment  of  1(X)  A.m.  The  fields 
for  this  source,  below  water  are  given  in  Bgure  5.  Again  the  fields  are  less  in 
the  three-layer  case  at  close  range  as  would  be  expected. 


5.  Conclusion 

The  three-layer  model  developed  herein  illustrates  the  importance  of  explicitly 
including  the  conductivity  of  the  sea-bed  in  calculations  of  the  magnetic  field 
from  a  horizontal  electric  dipole.  Specifically,  the  fields  could  be  up  to  50% 
less  than  might  otherwise  be  expected. 
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